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Abstract. Mathematical diffraction theory is concerned with the diffraction image of a given 
structure and the corresponding inverse problem of structure determination. In recent years, 
the understanding of systems with continuous and mixed spectra has improved considerably. 
Moreover, the phenomenon of homometry shows various unexpected new facets. Here, we report 
on some of the recent results in an exemplary and informal fashion. 



1. Introduction 

The diffraction measure is a characteristic quantity of a translation bounded measure a; on 
EucUdean space (or on any locally compact Abelian group). It emerges as the Fourier transform 
7 of the autocorrelation measure 7 of a;, and has important applications in crystallography, 
because it describes the outcome of kinematic diffraction (from X-rays or neutron scattering, 
say). In recent years, initiated by the discovery of quasicrystals (which are non-periodic but 
nevertheless show pure Bragg diffraction), a systematic study by many people has produced 
a reasonably satisfactory understanding of the class of measures with a pure point diffraction 
measure, meaning that 7 is a pure point measure, without any continuous component. 

Clearly, reality is more complicated than that, in the sense that real world structures will (and 
do) show lots of continuous components as well. Unfortunately, the methods around dynamical 
systems that are used to establish pure point spectra do not seem to extend to the treatment 
of systems with mixed spectrum, at least not in sufficient generality. Nevertheless, systems 
with continuous spectral components have recently been investigated also from a rigorous 
mathematical point of view, with a number of unexpected results. In particular, the phenomenon 
of homometry becomes more subtle, as we will see below. 

In this informal exposition, we summarise some classic results that have recently resurfaced in 
the context of diffraction, with special emphasis on singular continuous and absolutely continuous 
spectra. For proofs (and the formal details) we refer to the original papers or to work in progress. 



2. Mathematical setting and pure point spectra 

Below, we mainly consider various (weighted) Dirac combs on the real line, with support on Z 
(the set of integers), and their diffraction. Even in this simple and seemingly restricted setting, 
unexpected phenomena show up. If a; = with 5n the normalised point measure 

at n, the natural autocorrelation (if it exists) of u is 



(1) 



which we view as a measure on M. Here, = u;|[_7v,Ar] is the restriction of lo to the closed 
interval [— A^, A^], and p( is the 'flipped over' measure, defined by /1(g) = ^(g) with g a continuous 
test function of compact support and 'g{x) = g{—x). Here, and in analogous situations below, 
the bar denotes complex conjugation. We will only consider situations where the limit in ([1]) 
exists, either always (in the deterministic cases) or almost surely (in the probabilistic cases). An 
explicit calculation shows that supp(u;) C Z implies that 7 must be of the form 

with the autocorrelation coefficients 

1 ^ 

r?(m) = ^lirn^ 2A^ + 1 ^ w{n)w{n-m). (3) 

n=-N 

The existence of the limit in ([T]) is equivalent to the existence of the 77(771) for all ttt, E Z. 

By construction, 7 is a positive definite measure, and hence Fourier transformable. The 
result is 7, the diffraction measure, which is a positive, translation bounded measure on R. It 
describes (kinematic) diffraction from the measure compare [12j for background and jl8l [U] 
for the development of this approach to diffraction theory. One of the benefits of this approach 
is the unique decomposition 

7=(7)pp + (7)3c + (7)3c (4) 

of the diffraction measure into its pure point, singular continuous and absolutely continuous 
parts, the latter splitting relative to Lebesgue measure A, which is the natural reference measure 
for volume in Euclidean space. 

The best known example is the lattice Dirac comb u; = 5^ = ^^g^i^n, with autocorrelation 
7 = (5^ and diffraction 7 = (5^. Here, 7 follows from an elementary calculation, while the formula 
for 7 is a consequence of the Poisson summation formula 

6p = denser) 6p* (5) 

for an arbitrary lattice F (in M*^) of density dens(-r), with dual lattice F*, the latter defined 
by -T* := {x € M*^ | € Z for all y € F}; see [9] and references therein. The integer lattice is 
perhaps the simplest example for a system with pure point diffraction. More generally, if u; is 
any Z-periodic measure, it can be written as uj = fj, * 5^, with fi a finite measure. This leads to 
7 = {fi *]!) * 6^ and 7 = (5^, where is a continuous positive function on M. Its values at 
integer points are the intensities of the Bragg peaks. 

A similarly nice formula formula holds for regular model sets. Within a given cut and project 
scheme (CPS) (M'=',M™,/:) with lattice C C M''+'", consider A = {x € L \ x* € W}, where 
L = T^{C) is the projection of C in W^, W is the (compact) window in R™ (with boundary of 
measure 0), and -k denotes the star map of the CPS. By the general model set theorem 118 ^ 127 ^ 19]. 
the diffraction 7 of the Dirac comb is then the pure point measure 

^ =Y.\A{k)\'5u, (6) 

k£L® 

where L® = vr(>C*) is the projection of the dual lattice C* and the amplitude A{k) is given by 

dens(yl) [ orrik^y dens(yl) 



- vol(H^) - vol(H^) ^ > 



Table 1. Weights for the homometric pair from }17j. 
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This formula has various generalisations, for instance to model sets with other internal spaces 
or to weighted model sets; see \27\ [9] and references therein for more. 

As an example, we consider the period doubling sequence, as defined by the substitution 

g = ^jpd : a 1-^ ab , 6 i-^ aa. (7) 
A two-sided sequence can be obtained from a fixed point of via the iteration 

a\a abaa\abaa . . . — > w = q'^{w), (8) 

with convergence in the (obvious) product topology. Here and below, we write two-sided 
sequences as w = . . ■w_2W_i\wqWi . . . and use | to mark the origin. Note that a\a in ([8]) is 
a legal seed, so that w is a, fixed point of in the strict sense. 

We attach a Dirac comb to w hy uj = J2n£Z ^(^n) ^n, with h{a) = h+ and h{b) = This 
turns out to define a (weighted) regular model set [9] , so that the diffraction is indeed a pure 
point measure. It is given by 

7^= Y,\h+A{k) + h^B{k)f6„ 

where L® = lj£>i Z/2^ = {p | (r = 0, m G Z) or (r > 1, m odd)} is the Fourier module of the 
period doubling sequence. The amplitudes read 

Mk) = j-^e^^'' and B^k) = 5,^^ - A^k) , 

where we implicitly refer to the parametrisation of L®. We state this formula here without 
further details; see [ini E] for a proof. 

Let us briefly mention the homometry problem. It refers to the possibility that distinct 
measures can still possess the same autocorrelation. An interesting example is constructed in 
pT] . based on 6-periodic Dirac combs of the form * Yl^=o ^j^j- particular, the two choices 
of Table [Ulead to the same autocorrelation - and, in fact, even to identical correlation functions 
up to order 5. It is only the order 6 correlation that tells the two combs apart. 

For model sets, there are further possibilities, because different windows can have the same 
covariogram. The latter (for a compact set K) is defined as cvg^(x) = vol(i^ fl (x + K)^ = 

(l^ * l^)(x). This phenomenon results in identical autocorrelations and hence in homometric 
model sets. A simple planar example was constructed in ^; see also [1^I141[T3] . We skip further 
details here and shift our attention to continuous spectra now. 



3. Singular continuous spectra 

The paradigm of a singular continuous measure is the distribution function for the classic middle- 
thirds Cantor measure, which is also called the Devil's staircase. Here, we show a class of 
structures that lead to somewhat similar functions, yet with significant differences. 

The classic Thue-Morse (or Prouhet-Thue-Morse) sequence p!] can be defined via a fixed 
point of the substitution 

Q = £»tm: 1 i-^ 11 , 1 i-^ 11 (9) 
on the binary alphabet {1,1}. The one-sided fixed point starting with 1 reads 

= 1;q?;^|;2... = 11111111... , (10) 

while V is the fixed point starting with 1. One can now define a two-sided sequence w by 

]vi, i > 0, 
[v^i^i, i < 0. 

It is easy to check that w defines a 2-cycle under g, and hence a fixed point for g'^. Since the 
central seed 1|1 is legal, an iteration of g'^ applied to it converges to w in the product topology. 

The sequence w defines a dynamical system (under the action of the group Z). Its compact 
space is the (discrete) hull, obtained as the closure of the Z-orbit of w, 

^TM = {S^w \ i e Ij}, 

where S denotes the shift operator and the closure is taken in the local (or product) topology, 
where two sequences are close when they agree on a large segment around the origin. Now, 
(Xrp]y[,Z) is a strictly ergodic dynamical system (hence uniquely ergodic and minimal [25| [29]). 
Its unique invariant probability measure is given via the (absolute) frequencies of finite words (or 
patches) as the measures of the corresponding cylinder sets, which then generate the a-algebra. 

Here, we are interested in the diffraction of the (signed) Dirac comb Wrpj^ = ^nGz'^n^n, 
where we interpret 1 and I as weights 1 and —1. Its autocorrelation exists, as a consequence 
of unique ergodicity, and is actually the same measure for all sequences of Xrpj^. We can now 
employ the special structure of our fixed point w to analyse it. The autocorrelation is of the 
general form with coefficients rj(0) = 1, r]{—m) = rj{m) for all m € Z and 

r/(m) = lim — V VnVn+m (11) 
N^oo iV — ' 

71=0 

for all m > 0. Here, the structure of w and its relation to v was used to derive (jlip from ([3]). 
Observing that v satisfies = f„ and V2nJ^i = for all n > 0, one can derive from pT]) the 
recursions 

'q{2m) = 'q{m) and ry(2m -|- 1) = — — (r/(m) -|- r/(m -|- 1)) , (12) 

which actually hold for all m G Z. One finds ?7(±1) = — 1/3 from solving the recursion for m = 
and m = —1 with 77(0) = 1, while all other values are then recursively determined. 

To analyse the diffraction measure 7 of the TM sequence (following [23^ [20] ) , one can start 
with its pure point part. Defining S{N) = ^^_7v(^('^))^' ^'^^ derives via the recursion ()12p 
that S{AN) < f Z'(27V), which implies j^S{N) ) 0. By Wiener's criterion [30], this means 

(7) = 0, so that 7 is a continuous measure. 



Defining the (continuous) distribution function F via F{x) = 7([0, x]), another consequence 
of (I12p is the pair of functional relations 



dF(f)±dF(^) = {_J(^^^}dF(x). 



Splitting F into its sc and ac parts (which is unique) now implies backwards that the recursion 
(fT2|) holds separately for the two sets of autocorrelation coefficients, r]^^ and rj^^, with yet 
unknown initial conditions at 0. An application of the Riemann-Lebesgue lemma, however, 
forces ?/ac(0) = 0) hence rj^^{m) = for all m G Z, so that also (7)^^ = 0; compare [20j. 
This shows that 7, which is not the zero measure, is purely singular continuous. Figure [1] shows 
an image, where we have used the Volterra iteration 



1 f'^'' 

Fn+ii^) = 2/0 '^°^^'^y^)^^iy) with Fq{x) = X 

to calculate F with sufficient precision (note that F{x + 1) = F{x) + 1, so that a display on 
[0, 1] suffices). In contrast to the Devil's staircase, the TM function is strictly increasing, which 
means that there is no plateau (which would indicate a gap in the support of 7); see [5] and 
references therein for details. 

Despite the above result, the TM sequence is closely related to the period doubling sequence, 
via the (continuous) block map 

<f: 11, 11 a, 11, 11 6, (13) 

which defines an exact 2-to-l surjection from the hull X^pj^ to Xp^. 

The TM sequence is often considered as a rare and special example, which is misleading. To 
demonstrate the point (following [21]), consider the generalised Morse sequences defined by 

g = ek,i-- l-l'l', l-l'l' (14) 
for arbitrary A;,£ € N. Here, the one-sided fixed point starting with Vq = 1 satisfies 

_jv^, iiO <r<k, 



for m > 0. A two-sided sequence can be constructed as above. Since each choice of k, i 
leads to a strictly ergodic dynamical system, we know that all autocorrelation coefficients exist. 
Clearly (again with 1= — 1), we have 7]{0) = 1, while several possibilities exist to calculate 

k+£~3 

k+e+1- 



r7(ibl) = fjlxuf • In general, we obtain the recursion 



T]{{k + i)m + r) = {a^, g ., r?(m) + ak^^^k+i-r + 1)) , (16) 

with oij^i^ = k + £ — r — 2 min(/c, i,r,k + i — r), which holds for all m € Z and 0<r<k + i— 1. 

The recursion can once again be used to show the absence of pure point components (by 
Wiener's criterion) as well as that of absolutely continuous components (by the Riemann- 
Lebesgue lemma), thus establishing that each sequence in this family leads to a signed Dirac 
comb with purely singular continuous diffraction. The distribution function satisfies 



m>l 



F{x) =7([0,x]) =x+Yl 
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Figure 1. The distribution function of the TM measure on [0, 1]. 

which is a uniformly converging series. Moreover, the measure 7 has a (vaguely convergent) 
representation as an infinite Riesz product: With 6{x) := 1 + |r|^ X^^ii"^ '-'^fc £ r '^°^(^^^^)' 
is given by rin>o^((^ ~^ The entire analysis is thus completely analogous to that of the 

original TM sequence and shows that the latter is a typical example in an infinite family. Further 
details (and proofs) will be given in [3]. 

Let us finally mention that the block map (jl3p applies to any member of this family, and 
always gives a 2-to-l cover of the hull that belongs to the generalised period doubling substitution 

g': a^b^-\b^~^b, b ^ b''-^ab^~^a. 

Since we always have a coincidence in the sense of Dekking [13j, they all define systems with 
pure point spectrum - another analogy to the classic case k = i = 1. 

4. Absolutely continuous spectra 

The simplest example in this category is the Bernoulli (or coin tossing) comb. Let {Wn)n£Z be 
a family of independent and identically distributed (i.i.d.) random variables. We assume that 
W represents this family and takes values 1 and — 1 with probabilities p and 1 — p. We then 



consider the random measure 

^B = 5^W^n<5n, (17) 

which (almost surely in the probabilistic sense) has a natural autocorrelation of the form ([2]), 
with coefficients ^/^(O) = 1 and (a.s.) %(?ti) = {2p — 1)^ for all 7^ m G Z, so that we have 

7b = {2p - If 5^ + 4p(l - p) Jo (a.s.). (18) 

In particular, when p = 1/2, it gives 7g = 5^. Eq. (jlSp can be proved either by an application 
of Birkhoff's ergodic theorem for Z-action or by the strong law of large numbers (SLLN). We 
refer to [6] for details and further references. Let us also mention that (llSp can be rewritten 
^ 7b ~ /^^"^z "I" (1 ~ /5^)'^0' where /3 has the meaning of the drift velocity of a one-dimensional 
random walk on the line; compare [111 [8]. 

As a consequence, we obtain a diffraction measure of mixed type, 

% = {2p-lf5^ + Ap{l-p)X (a.s.), (19) 

where A denotes Lebesgue measure. The special case of a fair coin {p = 1/2) leads to Tb = A, 
which is purely absolutely continuous. More generally, when W takes the (possibly complex) 
values /i+ and each with probability 1/2, one obtains 

^ I hj-+h- |2 e- I h+—h- |2 , 

which is another simple example of a mixed spectrum with pure point and absolutely continuous 
components. 

An interesting deterministic counterpart is the Rudin-Shapiro sequence [T] . Its 

quaternary version is usually defined by the substitution 

Q = ^jpj^g : a 1-^ ac , 6 i-^' dc , c I— > a6 , d I— > db. 

A fixed point u = g'^{u) can be constructed from the legal seed b\a via iteration of g"^. The binary 
version follows then as the reduction w = f{u) with the mapping (p given by (p{a) = ip{c) = 1 
and 99(6) = ip{d) = —1. The two-sided sequence w satisfies 'w_i = —1 and Wq = 1 together with 
the recursion 

-(4- + ^) = b"^, . . for {0,1}, 
^ ' \(-l)"+^'u;(n), for £g {2,3}, ^ ^ 

for all n € Z. The sequence w and its hull once again define a strictly ergodic dynamical system, 
so that the natural autocorrelation 7pj^g exists and is the same for all members of the hull. It is 
of the general form To calculate it, one needs a little trick. Define the coefficients 



rj(m 



I := lim — w(n)w(n + m) I , 

n=—N 



(21) 



for m G Z. Note that all these limits exist, by Birkhoff's ergodic theorem (this is immediately 
clear for ri{m), but also follows for 'i?(m) on the basis of the quarternary sequence u, which also 
gives rise to a strictly ergodic system); compare |24] for an alternative approach. 

It is easy to check that 7/(0) = 1 and '!?(0) = 0, while a slightly more involved calculation 
reveals the following coupled system of recursions, 

1 + (-1)™ 

r/(4m) = r/(4m + 2) = 0, 

1 _ f-D™ f-D™ 1 

r/(4m + l) = ' r]{m) + ^ ^ i?(m) - ^ ^^(m + l), 

r/(4m + 3) = ' T?(m + 1) - ^ ^ ^9(m) + ^t?(m + l) , 



together with 

i?(4m) = 0, ??(4m+2) = ^ ^ + 1 7?(m + l), 

T9(4m + 1) = ^ r?(m) - ^ + ^ T?(m + 1), 

T?(4m + 3) = ^ r?(m + l)- ^ 7?(m) + ^ + 

These equations are vahd for all m € Z. Extracting = 1) = from the last two 
equations (with m = resp. ra = —1), one can then recursively conclude that T?(m) = for all 
m ^'L and r](m) = for all m ^ 0. This shows 

7rs = -^o and T^S = (22) 

In particular, the (deterministic) diffraction ()22|) of the binary RS sequence is the same as the 
(almost sure) diffraction of the coin tossing sequence ()19p for p = 1/2, despite the fact that their 
entropies are different (0 versus log(2)); see [121 E] for further details. 

The situation is actually 'worse' in the following sense. Let S G {±1}^ be a two-sided 
sequence, assumed ergodic, with associated Dirac comb a;^ = YlneZ'^nSn and autocorrelation 
7g. The latter exists due to the ergodicity assumption. Let {Wn)n<=z be the above family of i.i.d. 
random variables with values in {±1}. Consider now the new random Dirac comb 



Y^SnWndn, (23) 



which we call the BernouUisation of the original sequence. It can be considered as a 
'randomisation via second thoughts', since the action of the Wn is to either keep Sn (with 
probability p) or to change its sign (with probability 1 — p). By another (slightly more 
complicated) application of the SLLN, one finds [6] 

j = {2p-lf^s + M^-p)So (a.s.) (24) 

together with the corresponding almost sure diffraction 7 = {2p — 1)^7^ + 4p(l — p)X. If S is 
the binary RS sequence from above, one finds J = Sq and 7 = A, independently of the parameter 
p. This proves that the entire family defined by (j23|) for the RS sequence is homometric, with 
the entropy varying continuously between and log(2); compare [7] for the connection between 
diffraction and entropy in the pure point case. 



5. Outlook 

The diffraction measure is a useful tool, both for the understanding of experiments and for 
various theoretical questions. Beyond the well-studied case of pure point spectra, also continuous 
spectra are explicitly accessible and sometimes (as above) even computable. However, as was 
to be expected, various aspects are more involved, and this includes the inverse problem. The 
latter is increasing in complexity also with growing dimension. For instance, the Ledrappier 
system [22] in the plane has the same autocorrelation as the Bernoulli comb (both living on Z^), 
but has rank-1 entropy despite being genuinely two-dimensional. 

Various generalisations exist to non-lattice systems, which are best formulated via the theory 
of point processes p]. Here, the classic Poisson process (of mean point density 1) is another 
example (with random positions) that gives diffraction measure A when the points are randomly 
weighted with 1 and —1. It is clear that higher order correlations can tell these systems apart, 
but many questions are still open. 

Finally, more realistic models of randomness have to include interactions. The theory of 
Gibbs measures is a necessary and useful tool here, but only first steps in this direction have 
been taken; see [2] and references therein. Our present day understanding is still rather limited. 
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